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Abstract. In this article, we present a bifurcation analysis on the double- 
diffusive convection. Two pattern selections, rectangles and squares, are inves- 
tigated. It is proved that there are two different types of attractor bifurcations 
depending on the thermal and salinity Rayleigh numbers for each pattern. The 
analysis is based on a newly developed attractor bifurcation theory, together 
with eigen-analysis and the center manifold reductions. 



1. Introduction 

The main objective of this article is to develop a bifurcation and stability theory 
for the double-diffusive convection, using a newly developed bifurcation theory for 
dynamical systems, both finite and infinite dimensional [llj . Convective motions 
occur in a fluid when there are density variations present. Double-diffusive convec- 
tion is the name given to such convective motions when the density variations are 
caused by two different components which have different rates of diffusion. Double- 
diffusion was first originally discovered in the 1857 by Jevons [S], forgotten, and 
then rediscovered as an "oceanographic curiosity" a century later; see among others 
Stommel, Arons and Blanchard |T7], Veronis and Baines and Gill P]. 

We have conducted a bifurcation and stability analysis in 8] for the two di- 
mensional (2D) double-diffusive convection model. In this article, we continue our 
analysis to study the three dimensional (3D) double-diffusive convection problem. 
As in the two-dimensional case, the governing equations are the Boussinesq equa- 
tions with two diffusion equations of the temperature and salinity functions. 

The double-diffusive system involves four important nondimensional parame- 
ters: the thermal Rayleigh number A, the solute Rayleigh number 77, the Prandtl 
number a and the Lewis number r. We examine in this article different transi- 
tion/instability regimes defined by these parameters. We study the onset of double 
diffusive instabilities in these regimes. In the ocean circulation case, the Prandtl 
and Lewis numbers satisfy a > 1 and r < 1, as the heat diffuses about 100 times 
more rapidly than salt [16j . In this case, different regimes of stabilities and insta- 
bilities/transitions of the basic state can be described by regions in the X-ij plane 
(the thermal and salt Rayleigh numbers) as shown in Figure 13.11 In this article, 
we focus on the regimes where rj < rjc; see Figure IXTl In the case where rj > r]c, 
transitions to periodic or aperiodic solutions are expected, and will be addressed 
elsewhere. 
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In both 2D and 3D cases, the central gravity of the analysis is the reduction of 
the problem to the center manifold in the first unstable eigendirections, based on 
an approximation formula for the center manifold function. The key idea is to find 
the approximation of the reduction to certain order, leading to a "nondegenerate" 
system with higher order perturbations. The analysis for the 3D case studied 
here has all the difhculties appeared in the analysis for the 2D case such as the 
nonsymmetric linearized eigenvalue problem. Additional difhculties occur in the 3D 
case. One difficulty is the lack of the existence of global strong solutions of the 3D 
Boussinesq equations, and another is the higher dimensionality of the generalized 
eigenspace, leading to the center manifold reduction to higher dimensional center 
manifolds, and consequently to more complicated dynamics. 



With this reduction in our disposal, a general bifurcation theorem follows from 
the general strategy of attractor bifurcations. We prove in particular that there are 
two different types of transition/bifurcation from the basic state: one is subcritical 
and the other is supercritical. In the rectangle case (i.e. the horizontal domain is 
a rectangle), the types of transition/bifurcation are dictated explicitly by a nondi- 
mensional parameter r^^; see Figure [3?l1 In the square case (i.e. the horizontal 
domain is a square), types are determined by the quadratic form (j3.1ip . From the 
physical point of view, the supercritical bifurcation corresponds to a continuous 
transition. The subcritical bifurcation corresponds, however, to a jump transition, 
a very different transition. As demonstrated in the 2D case [S], this subcritical 
bifurcation leads to the existence of a saddle-node bifurcation and the hysteresis 
phenomena, thanks to the existence of strong solutions and the existence of global 
attractors. In the 3D case studied here in this article, we are not able to make 
such conclusions. However, we can still prove that there is an absorbing region 
resembling the main properties of the saddle-node bifurcation and hysteresis. 



There have been extensive studies on bifurcation and stability analysis for convec- 
tion models; see among others [T4j |2l IH Ej and the references therein. We mention 
in particular that in [3] , Buzano and Golubitsky used the singularity-theory method 
to study the problem of pattern formation as it relates to bifurcation with respect 
to the hexagonal lattice. They treated bifurcation problems that are symmetric 
with respect to group preserving doubly periodic functions on a hexagonal lattice. 
Their method provides a delicate alternative to study the convection problem in 
this paper. 

Another remark is that as mentioned by Sattinger on page 96 in [15| . a stable 
bifurcating solution with respect to a subclass of disturbances of full symmetry 
might turn out to be unstable with respect to the disturbances of full symmetry. 
In our analysis, although we use the symmetry properties implicitly to reduce the 
calculations, we don't use the equivariant theory to reduce the problem to a lower 
dimensional space. Since we want to keep as much stability information as possible. 

This article is organized as follows. The basic governing equations are given in 
Section 2. The main theorems are stated in Section 3. The remaining sections are 
devoted to the proof of the main theorems, with Section 4 on eigenvalue problems. 
Section 5 on center manifold reductions and the proofs of the theorems. 
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2. The Double-Diffusive Equations 



The nondimensional double-diffusive convection problem in a three-dimensional 
(3D) domain il = x (0, 1) C K."^ with coordinates denoted by {x, y, z) are given 
as follows; see Veronis 1191 : 



(2.1) 



^ = (j{AU - Vp) + (j{\T - T]S)e - {U ■ V){/, 
dT 

— = AT + w-{U-V)T, 
BS 

— = tAS + w~{U ■ V)5', 



div[/ ^ 0, 

where U = (m, v, w), X the thermal Rayleigh number, rj the salinity Rayleigh num- 
ber, a the Prandtl number, and r the Lewis number. We consider the periodic 
boundary condition in the x and y directions 

(2.2) (C/, r, S){x, y, z, t) = ([/, T, S)[x + 2j7r/a,,y, z, t) 

= ([/,T, S){x,y + 2kTT/a2,z,t), 

for any j, k € 1j. At the top and bottom boundaries, we impose the free-free 
boundary conditions; namely. 



(2.3) 

It's natural to put the constraint 



du 



{T,S,w)^0, — = 0, 



dz ' 



at z = 0, 1. 



(2.4) 



udxdydz 



vdxdydz = 



for the problem (|2.ip - (|2.3p . It is easy to see that (|2.ip is invariant under this 
constraint. The initial value conditions are given by 

(2.5) {U,T,S) ^ {U,f,S) at t = 0. 

Let 

H ={{U,T,S) e L^{n)^ \diyU = 0,w U=o,i= 0, (it,?;) satisfy ^ and (1^1. 
Hi =={([/, T, S) e ij2(rj)5 n I (f/, T, S*) satisfies (g^]) and 
Let G : iJi — > ff, and L^,, = — ^ — : Hi ^ H he defined by 

G(v) - (-P[(c/ • ^)u],-{u ■ v)r, -([/ • v)5), 

(-P[a(AC/)],-AT,-TA5), 
SAr,V' = (--P[CT(Ar - r/S')e], -w, -w), 

for any tp = {U,T, S) E Hi. Here P is the Leray projection to fields, and for 
a detailed account of the function spaces, see among many others [IB]. Then the 
Boussinesq equations (|2.ip - (l2.4p can be written in the following operator form 



(2.6) 



It is classical that there is a global weak solution for the system, and there is 
a global strong solution for small data; see among others [3 [TO]. Of course, the 
global existence of strong solutions is not known for large data. 
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3. Main Results 

3.1. Attractor bifurcation theory. In this subsection, we recapitulate the at- 
tractor bifurcation theory introduced by two of the authors in [T^l [H] . 

Let H and Hi be two Hilbert spaces, and Hi ^ H he & dense and compact 
inclusion. We consider the following nonlinear evolution equations 

(3.1) |^ = iv. + G(H,A), 

y m(o) = Uo, 

where u : [0,oo) H \s the unknown function, A G M is the system parameter, 
and L\ : Hi H are parameterized linear completely continuous fields depending 
continuously on A G M^, which satisfy 



(3.2) 



— Lx = A + B\ a sectorial operator, 
A : Hi H a linear homeomorphism, 
B\ : Hi —>■ H parameterized linear compact operators. 



It is easy to see [7 that L\ generates an analytic semi-group {e*^^}t>o. Then 
we can define fractional power operators (— L^)^ for any < M < 1 with domain 
Hfj, = D{{-LxY) such that H^,^ C H^j_^ if > //2, and Hq = H. 

Furthermore, we assume that the nonlinear terms G'(-,A) : H^ — > H for some 
1 > > are a family of parameterized C" bounded operators (r > 1) continuously 
depending on the parameter A G M^, such that 

(3.3) G{u,\)^o{\\u\\h^X VAgRI. 

In this paper, we are interested in the sectorial operator —L\ = A -\- B\ such 
that there exist an eigenvalue sequence {pk\ C and an eigenvector sequence 
{ek.hk} C Hi of A: 

'Azk=PkZk, Zk^Ck + ihk, 

(3.4) < Repfe ^ oo (A: ^ oo), 

^ \luipk/{a + Repfe)| < c, 

for some a, c > 0, such that {efe, /ife} is a basis of H . Also we assume that there is 
a constant < 6' < 1 such that 

(3.5) Bx -.He — >H bounded, V A G 



Under conditions (|3.4p and p.Sp . the operator —L\ = A+B\ is a sectorial operator. 

Let {S'A(t)}t>o be an operator semi-group generated by the equation (|3.ip . Then 
the solution of (|3.ip can be expressed as ^{t, ■0o) — S\{t)'iljo, t >0. 

Consider (|3.ip satisfying (|3.2p and (|3.3p . We start with the Principle of Exchange 
of Stabilities (PES). Let the eigenvalues (counting the multiplicity) of Lx be given 
by /?i(A), /32(A), • • • . Suppose that 

f < 0, if A < Ao 

(3.6) RePiiX) < = 0, if A = Ao {1 < i < m) 

[>0, if A > Ao 

(3.7) Rel3j{Xo) < 0. V m -I- 1 < j. 
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Let the eigenspace of L\ at Ao be 



oo 



En 



y \J{u,veHi\ {Lx„^ l3.j{Xo)fw^Q,w = u + iv}. 



l<j<m k—1 



It is known that dim_Bo = m. 

Theorem 3.1 (T. Ma and S. Wang [HIIII]). Assume that the conditions 

i3.5\) and i3.6\ )- (377^ hold true, and u = is locally asymptotically stable for iS.l]) 

at X — Xq. Then the following assertions hold true. 

(1) iS.l]} bifurcates from (u, A) = (0, Aq) to attractors S^, having the same 
homology as S""^^, for X > Xq, with m — 1 < dimTj\ < m, which is 
connected if m > 1; 

(2) For any Ux G u\ can be expressed as 



(3) There is an open set U d H with G f/ such that the attractor bifurcated 
from (0, Aq) attracts U\T in H , where T is the stable manifold ofu — with 
CO- dimension m. 

3.2. Main theorems. We now consider the double diffusive convection equations 
()2.ip . In this article, we always consider the case where the parameters A and rj 
satisfy 



Here the condition on ai and a2 defines the aspect ratios of the domain. The first 
eigenspace of the eigenvalue problem (|4.ip is of dimension two. 



Ux = vx + o{\\vx\\hi), vxeEo; 





Tic, Tic 



Figure 3.1. 



First we consider a more physically relevant diffusive regime where the thermal 
Prandtl number a is bigger than 1, and the Lewis number r is less than 1: 

(3.9) a>l>T. 
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In this case, we consider two straight hnes in the X — rj parameter plane as shown 
in Figure 13.11 

A = A,(,y) = ^ + 

T 4 

Also shown in Figure 13.11 are two critical values for rj given by 

r,^ = |7rV(l + a-i)(l - r)-\ = |.V(1 - r^' ■ 

The following two theorems study the transitions/bifurcation of the double- 
diffusive model near the line Li for t] < rjc- 

Theorem 3.2. Assume that the condition i3.8\) holds true, a > 1 > t, and rj < 
rjc^ — ^t:^t^{1 — T^)^^. Then the following assertions for the problem i2.1\) - [2^\ l 
hold true. 

(1) If X < Ac, the steady state {U,T,S) = is locally asymptotically stable for 
the problem. 

(2) For X > Xc, the solutions bifurcate from {{U,T, S), X) = (0, Ac) to an at- 
tractor ~ , consisting of steady state solutions of the problem. 





H 








V ''' 



X 



Figure 3.2. If t > 1 or 77 < 77C1, the solutions bifurcate from 
(0, Ac) to an attractor Sa for A > Ac. 



Theorem 3.3. Assume that the condition LS. 8\) holds true, a > 1 > t and rjc > rj > 
rjc^ — ^n^T^{l — T^)^^ . Then on X < Xc, the problem i2.1\) - [K^ bifurcates from 
(([/, T, 5*), A) = (0, Ac) to a repeller E>, = , consisting of steady state solutions of 
the problem. 
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H 








X ] /x^ 



Figure 3.3. If r < 1 and rjc^ < i] < rjc, the problem has a sub- 
critical type bifurcation for A < Ac. 

We now consider the diffusive parameter regime where cr > 1, r > 1, and a =/= t. 
In this case, two lines are shown in Figure 13.41 The following theorem provides 
bifurcation when A crosses the line Li. 

Theorem 3.4. Assume that a>\,T>\,a^T and \3.^] hold true, then for any 
rj > 0, the following assertions for the problem \2.1\) - (2^ hold true. 

(1) If X < Ac, the steady state {U,T,S) = is locally asymptotically stable for 
the problem. 

(2) For A > Ac, the solutions bifurcate from {{U,T, S), X) = (0, Ac) to an at- 
tractor Sa, homeomorphic to , which consists of steady state solutions of 
the problem. 




n 



Figure 3.4. 

3.2.2. Square Case. For the square case, we assume ai = a2 = f • Then, the 
first eigenspace of the eigenvalue problem (|4.ip is of dimension 4. In this case, the 
bifurcation type is determined by the following quadratic form 

(3.11) (5i(a;iiii + + a;?_ni + yl-mf + 2(52(a;nn + ynii)(2:?_ni + 

where Si and 82 are as defined in (|5.7p and (|5.8p . which can be evaluated precisely 
when the parameters r, cr, rj and A are given. More precisely, we have the following 
two theorems. 
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Theorem 3.5. Assume that i3.8\) holds true and liS.ll]) is negative definite, then 
the following assertions for the problems li2.1\} - (K4\ l hold true. 

(1) If X < Ac, the steady state (U, T, S) = is locally asymptotically stable for 
the problem. 

(2) For A > Ac, the solutions bifurcate from {{U,T, S), X) — (0, Ac) to an at- 
tractor having the same homology as S'^ . 

Theorem 3.6. Assume that i3.8\) holds true and iS.ll]) is positive definite. Then 
on X < Ac, the problem i2.1\) - (K^ bifurcates from {(U,T, S), X) — (0, Ac) to a 
repeller E^, having the same homology as S^. 



4. Eigenvalue problem 



In order to apply the center manifold theory to reduce the bifurcation problems, 
we consider the following eigenvalue problem for the linearized equations of (|2.ip - 
(Ell): 

a{AU - Vp) + a{XT - i^S)e = PU, 
AT + w = f3T, 
tAS + w = /3S, 
divU = 0, 



(4.1) 



supplemented with (|2.3p and 



4.1. Eigenvalues. We shall use the method of separation of variables to deal with 
the problem (|4.ip . Since i/j = {U,T,S) satisfies the periodic condition (|2.2p . we 
expand the fields in Fourier series as 



(4.2) 



j.k— — oo 



(4.3) 



Plugging (|4.2p into (|4.ip . we obtain the following ODE system: 

DjkUjk - ijaipjk = a^^Pujk, 
DjkVjk - ika2Pjk = cr^^Pvjk, 
DjkWjk ~ p'jk + XTjk - r]Sjk = a^^Pwjk, 
DjkTjk + Wjk — PTjk, 
rDjkSjk + Wjk = f3Sjk, 
ijaiUjk + ika2Vjk + w'^f. = 0, 
^ w^-fe U=o,i= v'jk U=o,i= Wjk \z=aA— Tjk \z=o,i~ Sjk \z=o,i~ 0, 



for j,k e Z, where ' = d/dz, Djk — d^/dz^ — a^^, and a^^, — j^a\ + fc^a^ 
Wjk 7^ 0, (|4.3p can be reduced to a single equation for Wjk{z): 

(4.4) {{TDjk - l3){Djk - P){Djk - <y-^P)Djk 



If 



(4.5) 



Wjk 



+ a%[X{TDjk - /?) - viDjk - P)]}wjk = 0, 



Wjk = w 



jk 



= wjl' =0 at z = 0, 1, 
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for j, k Thanks to (|4.5p . Wjk can be expanded in a Fourier sine series 

oo 

(4.6) Wjk {z) ^ Wjki sin Inz, 

1=1 

for {j,k) e Z X Z. Substituting (|4.6p into (|4.4p . we see that the corresponding 
eigenvalues /3 of Problem (|4.1[) satisfy the cubic equations 

(4.7) f3' + {a + T+ l)7|fez/32 + [{a + r + aT)^%, - aa^lji^^iX ~ ry)]/? 

+ + cr"jfe('7 - tA) = 0, 

for j,k gZ and Z G N, where 7^^,; = + ^^tt^. 

For the sake of convenience to analyze the distribution of the eigenvalues, we 
now introduce some notations. For fixed parameters a, t, rj and A, let 

gjki (13) = + (a + T + l)7'fe//3' + ((T + T + ar)7|fe,/3 + (7t-/%i , 

hjki{l3) = [fTa,\77fe7(A - - <^(^%iv - rX), 

27 

^c.=f^V3(l_rr\ A. = ^ + |7r^ 
Furthermore, let (3jkii, Pjki2 and /3jfei3 be the zeros of fjki with 
ReiPjui) > ReiPjkn) > ReiPjkis), 

and let (3jkio- = 

The following lemma characterizes the eigenvectors of L\ri with zero w-components; 
the proof is straightforward and we omit the details. 

Lemma 4.1. For / G U N and {j, fc) G Z^, Pjkia is an eigenvalue of the problem 
^■1^ . Moreover, the following assertions are true: 

(1) // (j, k) — (0, 0) and I > 0, the corresponding eigenvectors are 

(4.8) Vf™" = (cos It:z, 0, 0, 0, 0)*, i/'2™"' = (0, cos Zttz, 0, 0, 0)*; 

(2) If + fc^ 7^ and ^ G U N, the corresponding eigenvectors are 

ipf^"" — (k cos{jaix + ka2y) cos Zttz, — j cos(jQ!ia; + ka2y) cos Zttz, 0, 0, 0)*, 

(4.9) 0. . 

^^3kic7 _ s\n{jaix + ka2y) cos ^ttz, — j sin(jaix + ka2y) cos Zttz, 0, 0, 0)*. 

In the following discussions, we shall focus on the following diffusive regime: 

(4.10) (T>l>r>0, r] < rjc and A w Ac. 

Lemma 4.2. Suppose that a^j, — ^ and (^T7^ /loWs true, then fjkiiP) has three 
simple real zeros. 

Proof. Since A w Ac, it suffices to prove this statement for A = Ac. In this case, 
IMP) + [37r'('^ + r + l)/2]/3' + [QnHa + r + ar)/4 - a{X, - ry)/3]/3 
=/(/3)/3. 

Since rj < rjc, the constant term of /(/3) is positive. Hence f3jkii = is a simple 
zero of fjki- Moreover, the quadratic discriminant of /(/3) is 97r"'((T + 1 — t)^/4 + 
Aar]{l — t)/{3t) > . This implies fjki has three simple real zeros. □ 
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We summarize the following important lemma about the distribution of the 
zeros of fjki- From the physical point of view, this lemma verifies the principle of 
exchange of stabilities (PES). 

Lemma 4.3. Assume that either 1) < rj < rjc with t < 1 or 2) rj > with 
T > 1 holds true, then 

{<0 if X< Xc 2 
= */ A = A, z/ a% = ^, 

> if X> Xc 

(4.12) Re(3,ki,iX) < */ {a%,l, q) ^ (^, 1, 1). 

Before proceeding to the proof, we make a few remarks as follows. 

Remark 4.4. The distribution of the zeros of fjki was first partially analyzed by 
Veronis [T^; see also P.G. Baines [T]. The complete proof of this lemma is similar 
to the 2D case in [S], we shall omit it. 

To check that the operators —L\n satisfy condition p.4p . we prove the following 
lemma. 

Lemma 4.5. (1) All but finitely many zeros of fjki{P) are negative real num- 
bers for (j, k, Z) e X N. 

(2) /3jkiq -oo iff + k^ + P -^(X. 

Proof. Since fjki = gjki — hjki, /? is a zero of fjki{(3) if and only if f3 satisfies the 
equation 

(4.13) gMl3) = hMl3)- 
Plugging P = "f^^^iP* into KW . we obtain 

(4.14) {13* + l)(r + t)(/3* +a)= d.kdiX ~ - iv ~ rX)], 

where 'djki = <^<fjkll%i- Since limp^k^^ii^^'djki = 0, the roots of (|4.14p must be 
negative real numbers near the interval [— cr, — r] when {f + fc^ + P) is large. This 
completes the proof. □ 

4.2. Eigenvectors. We now make some observations to analyze the spectrum of 

Lxrr Since 5jfc/(/?) = (^ + 7|fe;)(/3 + T-i%i){l3 + a-f^i) and hju = (^a^lJkAi^ - 
v)P ~ iv ~ '''^)ljki]' it's 6^sy to check that (3 = —Jjki or (3 — —T^'jki is a zero of 
fjki[f3) if and only if ajk = 0. In the case of a^j, — 0, the zeros of /oo; are — t7qq; 
(/3oo/i), -7ooz (/3oo/2) and -cr7oo/ {Pmz)- The corresponding eigenvectors are 

V''^™"(a;,j/,z) = (0,0,0,0,sinZ7rz)*, 1^^°°'^ {x , y , z) = (0, 0, 0, sin^Trz, 0)*, 

(4-15) 

V^f {x, y, z) = (cos Zttz, 0,0,0,0)*, V2 " i^^ ^) = (0, cos ^tt^, 0, 0, 0)*. 
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To analyze the structure of the eigenspaces of (|4.ip . for fc G Z, j e {0} U N and 
Z e N, we define 

9-iki — \ — 9 — cos( jaia; + ka2y) cos Lttz, — ^ — cos( jaix + ka2y) cos jttz, 

sinQ'aia: + ka2y) sin^Trz, 0, 0)*, 
4''jki = (0: O7 0, sin(jQ;ia; + ka2y) sin Zttz, 0)*, 
'/'jfei = (0, 0, 0, 0, sin(jQ;ix + fcQ;2y) sin Zttz)*, 

4 jailn ... , ka2lTT , . . , 
= ( 2 — sin(jaix + ka2y) cos «7rz, ^ — sin(jQ;ia; + ka2y) cos /ttz, 

cos(jq;ix + ka2y) sin /ttz, 0,0)*, 
— (0, 0, 0, cos(jQ;ia; + fcQ;2?/) sin/TTZ, 0)*, 
•^jfci = (0' 0' 0' 0, cos(jQ!ix + ka2y) sin Zttz)*. 
The following lemma follows from (|4.ip - (l4.7p . 

Lemma 4.6. If + fc^ 7^ anrf PjUq (q — l,2,3j is a zero of fjki, then we have 
the fallowings. 

(1) The eigenvector corresponding to (3jkiq in the complexified space of H is 
(4.16) i;"^'"^ ^e^(J"^"+'="-^^)(^:^^cosZ^z,^^^cos?7rz,sin/^z, 

^1 (/3j7cig) sin Zttz, A2{/3jkiq) sin Zttz)*, 

(2) // is a real number, the corresponding eigenvectors are given by 

^2^"' = + ^1 {M4>%1 + A2 {Pjklq)c^%i 

(3) If Im{Pjkiq) 7^ 0, the generalized eigenvectors corresponding to (3jkiq and 
Pjklq are 

(4.18) 

^1"'' = <P]kl + Rl{Pjklq)cf]kl + R2{(3jklq)<P%i + h{Pjklq)4>%i + HP ,klq)4>%i . 
V,f^"' = -/l(/3jM,)<^,\z - l2{fijklq)4>]kl + <t>%i + i?l(ftM,)<^ffeZ + R2{Pjklq)<l>%u 

4""' = 4>]ki + Ri{hkiqWjki + i?2(,3,M,)0,\i + h{hkiq)<i>%i + /2(^jfe;9)</>,\z, 

V's^^"" - -imklq)^]kl - Hkklq)4>%i + + Rl{kklq)<l>%i + (^, feig ) M , 

w/iere i?i(/3) = i?e(Ai(/3)), /i(/3) = Im{Ai{l3)), i?2(/?) = -Re(A2(/3)) and 
/2(/3) =/m(A2 (/?)). 

Definition 4.7. (1) //j = fc = 0, for each I eN, we define Eqqi = spa7i{?A^oo!i ^ ^/3ooi2} 

and E^Qi = spanlV-f™'', '02™"'}- 
(2) Forf+k"^ ^ 0, we define E]^^i = spa^^jfci, '/'jfez, '/'jfe;}: ^^|fci = •spa«{<^|fc;, 0^^;}, 

Ejki = ii^jfci © and i?;,, = WV^f^"'-, 
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(3) For (j, fc, Z) G X N, we define Ef.^^ to be the eigenspace spanned by the 
eigenvectors and the generalized eigenvectors corresponding to the zeros of 
fjkl ■ 

Remark 4.8. (1) It is easy to see from the Fourier expansion that {Ejki U 
Ejki}Ti=o,k=-oo is a basis of Hi. 
(2) Ejki (resp., EJ^^^) is orthogonal to E^.k^h ("^(^^P-^EJ^k^h) ^r {j,k,l) ^ 
(ji, ki,li), and Ejki is always orthogonal to EJ^f.^1^. 

The following theorem together with Lemmas 14.11 and 14.31 complete the analysis 
of the eigenvalue problem (|4.ip . 

Theorem 4.9. Under the assumption ^.10^ , we have 

1) Ef^^, = Ejki for j e {0} UN, ke Z, and I £ N; and 

2) Lxf_i\Ejki is strictly negative definite for each {j, /c, Z) G x N when A < Ac- 

Proof. The first assertion follows by the fact that Ejki is an invariant subspace of 
Lx,i and fjki is the characteristic polynomial of La?; \Ejki ■ Assertion 1), Lemma 14.31 
together with the fact that fjki does not have a zero of multiplicity three imply 
Assertion 2). 

□ 



1 oc 

lj,l=0,k= 



We conclude the above analysis as follows.. 

(1) The eigenvalues of Lxrj : Hi ^ H consist of {Pjkiq, PjkiaVj 

(2) The (generalized) eigenvectors of Lx,^ form a basis of H. 

(3) —Lxri is a sectorial operator. 

(4) For the rectangle case, the multiplicity of the first eigenvalue, /3ioii('^), is 
2, and the corresponding eigenvectors are and 

(5) For the square case, the multiplicity of the first eigenvalue, /3iiii(A)(=/3i_iii(A)), 
is 4, and the corresponding eigenvectors are V'f^'^'^S ^2^^^^, t/j^^""^" and 

To show that G satisfies condition p.3p . We assume | < /i < 1, then for 
e H^ C H, by Sobolev's inequality, 

\Gmi < f r'"^ r'"' \M^m\'dxdydz < mumi < c\^p\%^. 

Jo Jo Jo 
where C is some constant. Hence, G{tp) = o(|-!/)|/f^). 

4.3. Dual Basis. Since Ejki is finite dimensional for each {j,k,l), there exists a 
vector ^p^*" e Ejki [q = 1, 2, 3andp = 1, 2) such that 

= if {q*,p*) ^ {q,p). 



We choose 



Hence, by the orthogonality of Ejki and EJ,^i, {*p''", *p'""}°^,=o,fc=-oo fo™ a 
dual basis of H corresponding to {Ejki, EJi^iJ^^i^^ k=-oo i'^ the sense that 

=0 if (/3*,p*)^(/5,p). 
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The proof of Lemma 14.51 implies that aU but finitely many of fjki have three 
distinct real zeros. For such fju with + P ^ 0, vp^^'"''' and ^i^'''''' could be 
chosen as 



where 
(4.21) 



5. The Proofs of the Main Theorems 

We are now in a position to reduce equations of (|2.ip - (|2.4p to the center manifold. 
We would like to fix 77 < rjc, and let A « Ac be the bifurcation parameter. For any 

il) = {U, T, S) e H, we have 

00 00 3 00 00 

^= E E E(^j'='W^^"+2/.hX^"')+ E E(^j'«'-v^i'""+%'='-^2'"" 

k— — ocj—0,l—lq—l k— — ooj,l—Q 

For the square case, we assume ai = 0:2 = f . Hence /Sim (A) = /3i_iii(A) are the 
first eigenvalues. 

For brevity, we let /3o(A) — /3iiii(A) = /3i_iii(A), then the reduced equations 
are given by 



(5.1) 



— - /?o(A).i-ni + ^ ^ < Gi^,y^),^m >H, 



dyi-ni . , 1 



^ = /3o(A),i_ni + < G{^,^),^^m >H . 

Here for ^Pi = (f/i,Ti,^i), i/;2 = {U2,T2,S2) and ^3 = (C/3, 73,^3), 

pi p27v/a2 p27v/ai 
<G(Vl,V2),V'3 >if=- / / / [< (C/i • V)C/2,C/3 >M3 

Jo Jo Jo 
+ (Ui ■ V)T2Tz + [Ui ■ V)S2S3]dxdydz. 

Let the center manifold function be denoted by 
(5.2) 

$ = E (^'i(a;iiii, 2/1111, a;i_iii,yi_iii)-0f +$^(xiii,?/iii, a;i_iii,yi_iii)V'2 )■ 

/3#/3iiii,/3i_iii 

Note that for any G Hi (i = 1, 2, 3), < G('0i, i/'2), V'2 >h= 0, < G(V'i, V'2), V's >h= 
- < G{i>i,ip3),ip2 >H, and for any e E'jfe/ (i = 1,2,3), < G(V'i, V'2), i/'3 >i/= 0. 
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Applying Theorem 3.8 in [TT], we obtain 

$/3oo.i = + yl,,^ + a^tiii + t/tiii] + o{2), 

$/3oo.. ^ + + xtiii + Vl-iii] + o{2), 

= B(/32029)(a;iiii2/i-iii + a;i-iiiyiiii) + o(2), 
$f = B(/3202,)(-a;iiiia;i-iii + t/inij/i-in) + o(2), 
= i?(/3o22g)(a;iiii2;i-iii - a;i-iii2/iiii) + o(2), 
= B(/3o22g)(a;iiiixi„in + yiiiiyi-m) + o(2), 

(5.3) $(2^1111, J/llll,Xl_ni,yi-lll) = $/30021^/30021 +$/30022^/30022 

3 

+ ^(^f °''Vf'°'' + ^-f'^'^Vf""' + ^^'''"'Vf''^" + $2"'''V^2"''') + 0(2), 
9=1 



where 



S(/3) 



^(3 + Ai(/?o)Ci(/3) + A2(/3o)C2(/3)) 
' 2/3(5 + Ai(/3)Ci(/3) + A2{l3)C2m ' 



o(2) =o(2;?iii + yJiii + + yi^ni) 

+ 0(|/3o(A)| • (a:Lii + ynii + ^^Liii + yLiii))- 
Hereafter, we make the fohowing convention, 

o(3) =o((x2iii + yf^^^ + x?_iil + yf-iuf^^) 

-111 + J/i-iii 

)3/2)^ 

o(4) =o((a;?iii + y^i^ + + y?_ni)^) 

+ 0(|/?o(A)| • (a;?,!! + y^^^ + a^^iii + yLiii)')- 
By ()5.3p and the fact that /3o22(} = P202q^ we obtain 

< G(^,V),^'^'"" >H= < G(^f"",$),^^'"" >H + < G(V'^"\$),^^'"" >ff 

+ < G(<^-"\$),*f"" >H + < G(V'^^-"^),$),*f"" >H +o(3) 
= - < G(7/;f"",^'^'""),$ >_ff - < G(i/.f"",«'^'""),$ >ff 
- < G(<-^",*f""),$ >H - < G(V'2^^-"^),vI/f""),ci> >H +o(3) 

2,AiCi A2G2-, 2,2,2 , 2 N 

=271- (-5 + -5 )(a;iiii + yiiii + + yi_in)xiiii 

P0022 P0021 

3 

+ 2^ B*(/3202,)B(/3202,)(a;i-iii + y?-iii)a;iiii + o(3), 
9=1 

where 

(5.5) B*{p) = n{3 + Ai{p)Ci{Po) + A2(/3)G2(/3o)). 
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Similarly, we derive that 

< Gii;, ^t'' >H^27:\^ + ^)(x?ni + ^nii + ^tiii + yLiii)2/iiii 

P0022 P0021 

3 

+ 2^ B*(/3202g)B(/32029)(x?„ni + 2/Liii)yiiii + o(3), 
9=1 

< G(V, V), vl'f-- >H=2n^i^ + + yl,,, + xtni + 2/?-iii)^i-iii 

P0022 P0021 
3 

+ 2^ B*(/3202g)-B(/32029)(x?iii + ^I'lii )a:i-iii + 0(3), 

9=1 

< G(^, vl^f-"^ >^=27r2(4i^ + ^){x\,,, + y^^ii + xtni + 2/?-iii)yi-iii 

P0022 P0021 
3 

+ 2^ B*(/32029)5(/3202,)(x?ni + y?iii)2/i-iii + o(3). 

9=1 

Therefore, (|5.ip can be rewritten as 



(5.6) < 



rfs^iiii ^ /3(,(A)xiiii + (5i£;a;iiii + 52{x\_^^-^ + ?/i_iii)a;iiii + o(3). 



dt 



-^-^p- = /3o(A)yiiii + SiEyiiii + 52(a;?_iii + 2/?_iii)yiiii + o(3), 

'^^^^""^ = /3o(A)xi_iii + 5iExi-iii + (J2(a;nii + 2/iiii)a;i_iii + o(3), 

= /3o(A)yi-iii + (5i£;yi-iii + 52(a;nii + 2/nii)yi-iii + o(3), 
where E = (.t^^ii + y^^^ + x^^i^ + 2/?-iii), 

(5.7) 5i =Ji(A,ry) =2^2(^_^^)/4(3^^^^^^^^^^) 

P0022 P0021 
(AiCi +T-IA2C2) 



8(3 + ^1^1+^2^2)' 
3 

(5.8) (^2 =<52(A,77) -2^S*(/3202g)S(/3202g)/4(3 + AlCi+A2C2). 

9=1 

The energy estimate of (|5.6p is given by 
1 ^iii/ 

(5.9) -— - /3o(A)i? + S^E^ + 2J2(a;?ni + y?iii)(2:?-iii + ytiii) + "(4). 
If the quadratic form 

(5.10) 5ii;2 + 2,52(a;iiii+y?iii)(a;? -Ill + 2/i-iii)j 

is negative definite at A = Ac, we conclude that {U,T, S*) = is a locally asymp- 



totically stable equilibrium point of (|2.ip - (|2.4p . We then obtain Theorem 13.51 bv 
Theorem 13.11 directly. And if the quadratic form (|5.10p is positive definite, we 
obtain the subcritical bifurcation. This completes the proofs of Theorem 13.51 and 
Theorem [ 
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For the rectangle case, the first eigenspace is of dimension two. The reduced 
equations on the center manifold are given by 

dxiou 



(5.11) 



= /3ioii(A)xioii + (^i(a;ioii + yioii)a;ioii + o(3), 



dt 

dywii _ p^^^^(^x^y^^^^ _,_ Siixlf)^-^^ + yloii)ywii + o(3). 



dt 

The rest part of the proofs of Theorem 13. 21 Theorem l3.3l and Theorem 13.41 is similar 
to the 2D case in [8]. In this case, we can use Theorem 5.10 in [13] or invariant 
sphere theorem in [4] to conclude the structure of the bifurcating solution. 



Remark 5.1. (1) As we have already proved in [8], we know that 
(5.12) 6i 



< if 7] < , 

> if 77 > ?7ci ■ 

(2) When cr, r, A and 77 are given, 5i and S2 can be evaluated numerically. 
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